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Abstract 

A physical explanation for the saturation of broadband shock-associated 
noise (BBSAN) intensity with increasing jet stagnation temperature has 
eluded investigators. An explanation is proposed for this phenomenon with 
the use of an acoustic analogy. To isolate the relevant physics, the scaling 
of BBSAN peak intensity level at the sideline observer location is exam- 
ined. The equivalent source within the framework of an acoustic analogy for 
BBSAN is based on local held quantities at shock wave shear layer inter- 
actions. The equivalent source combined with accurate calculations of the 
propagation of sound through the jet shear layer, using an adjoint vector 
Green’s function solver of the linearized Euler equations, allows for predic- 
tions that retain the scaling with respect to stagnation pressure and allows for 
saturation of BBSAN with increasing stagnation temperature. The sources 
and vector Green’s function have arguments involving the steady Reynolds- 
Averaged Navier-Stokes solution of the jet. It is proposed that saturation 
of BBSAN with increasing jet temperature occurs due to a balance between 
the amplification of the sound propagation through the shear layer and the 
source term scaling. 
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i Personal Introduction 


2 It is a privilege to contribute an article to this special edition in honor 

3 of Dr. Fereidoun ‘Feri’ Farassat. The present article involves the use of 

4 an acoustic analogy and a Green’s function for its solution. Dr. Farassat ’s 

5 career was heavily involved with both of these fundamental methods in aeroa- 
e coustics since his Ph.D. [1] work at Cornell (under advisement of Professor 

7 William R. Sears) based on the work of J. E. Ffowcs Williams and D. L. 

8 Hawkings [2], His Ph.D. work laid the foundation for the rest of his career at 

9 NASA Langley Research Center (LaRC) within the Aeroacoustics Branch. 

10 Dr. Farassat ’s developments such as Formulation 1 (Farassat [3]), Formula- 

11 tion 1A (Farassat and Sued [4]), the use of generalized functions (Farassat 

12 and Myers [5]), and countless others, were extremely important for the field. 

13 Some of these contributions are available on the NASA Technical Reports 

14 Server, where Dr. Farassat has over 130 publications available to the public 

15 on a wide range of topics. 

io Dr. Farassat, during his mid- to late-career, was undoubtably the the- 
17 oretical backbone of the Aeroacoustics Branch at NASA Langley. He had 
is influenced the technical direction of many researchers within both the branch 

19 and NASA as a whole, and had a considerable influence throughout the com- 

20 munity, all of which are still being felt today. 

21 Dr. Farassat had a long history of imparting his knowledge to new re- 

22 searchers at NASA Langley. Some of my first and most memorable interac- 

23 tions with Dr. Farassat had started with these teachings. I enjoyed many 

24 technical discussions in his office and his guidance changed my technical view- 

25 point, especially relating to the acoustic analogy. These discussions saved me 

26 large amounts of time and helped me avoid many possible technical failures. 

27 He also was not afraid to offer advice, technical or personal, and was gen- 

28 uinely interested in the well-being of everyone he interacted with. He was an 

29 unwavering advocate within NASA for the importance of research and was 

30 extremely supportive of junior researchers. 

31 I am proud to call Dr. Farassat my colleague and friend. Thank you Feri 

32 for the time we had together. 

33 1 . Introduction 

34 Unfortunately, there is no first principles mathematical model or physi- 

35 cal understanding of how broadband shock-associated noise (BBS AN) scales 
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36 with increasing stagnation temperature. This paper attempts to examine the 

37 scaling of BBSAN intensity with increasing stagnation temperature via the 

38 acoustic analogy of Morris and Miller [6]. This is accomplished by examining 

39 the peak intensity at the sideline location relative to the jet centerline axis. 

40 The equivalent source of the BBSAN is modeled with the use of local instead 

41 of ambient quantities of a steady Reynolds-Averaged Navier-Stokes (RANS) 

42 solution of the jet exhaust and a simple model of the two-point velocity 

43 cross-correlation. Noise propagation is accurately modeled by using an ad- 

44 joint vector Green’s function solver for the linearized Euler equations (LEE). 

45 The scaling is compared with the measurements of Kuo et al. [7] for a design 

46 Mach number = 1.50 nozzle at over- and under-expanded conditions and 

47 with the measurements of Bridges and Brown [8] for a convergent nozzle. 

48 Comparisons cover the range of total temperature ratios ( TTR ) from one 

49 to four. The equivalent source model combined with accurate calculations 

50 of the propagation of BBSAN through the jet shear layer allows for predic- 

51 tions that retain the scaling with respect to nozzle pressure ratio (N PR) and 

52 allows for the saturation of BBSAN with increasing TT R. 

53 Jet noise is due to multiple unique sources. Far-held lossless noise spectra 

54 from an off-design singlestream supersonic jet can be observed in the far- 

55 held as shown in Fig. 1. The x-axis represents non-dimensional frequency as 

56 Strouhal number, St, which is frequency normalized by the fully expanded 

57 jet velocity, iq , and the fully expanded jet diameter, Dj. The y-axis repre- 
ss sents the Sound Pressure Level (SPL) per unit St referenced to twenty micro 

59 Pascals. The observer angle tfj is measured from the upstream axis of the jet 

60 centerline to the observer in the far-held about the nozzle exit plane. The 

61 non-dimensional distance from the nozzle exit to the observer is R/D = 100, 

62 where R is the distance and D is the nozzle diameter. 

63 Shock-associated noise consists of discrete tones often called ‘screech,’ hrst 

64 observed and described by Powell [9] ; and BBSAN was hrst extensively mod- 
es eled and studied by Harper-Bourne and Fisher [10]. Screech (see Raman [11] 
ee for an overview) has a large effect on BBSAN which will be illustrated later. 
67 BBSAN results when large-scale coherent turbulent structures interact with 
es the shock waves in the jet shear layer. Each interaction of turbulence with in- 

69 dividual oblique shock waves represents a source that contributes to BBSAN. 

70 The noise combines constructively or destructively in the far-held to produce 

71 the broad humps that are seen in Fig. 1. BBSAN is less intense in the down- 

72 stream direction than mixing noise due to refraction effects. In the sideline 

73 and upstream directions (i[) = 90 and ^ = 50 deg. in Fig. 1 respectively) BB- 
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74 SAN dominates the mixing noise over a wide range of frequencies. The peak 

75 frequency of BBSAN varies with observer angle and jet operating conditions. 

76 For overviews of jet noise consult Ffowcs Williams [ 12 ], Ffowcs Williams [ 13 ], 

77 or Goldstein [ 14 ]; and specifically for supersonic jet noise consult Tam [ 15 ]. 

78 The question arises regarding how BBSAN scales with increasing tem- 

79 perature. The scaling of BBSAN with increasing TTR can be observed 
so experimentally in Fig. 2 . The trend is similar (in terms of intensity scal- 
si ing and ‘saturation’) across observer angles, jet Mach numbers, and nozzle 

82 geometries. The phenomenon is summarized excellently by Viswanathan et 

83 al. [ 16 ] who state, “The levels increase as the jet is first heated; however, the 

84 levels do not increase with further increase in jet temperature. The physical 

85 phenomenon responsible for this saturation of levels is not known at this 

86 time.” 

87 Harper-Bourne and Fisher [ 10 ] observed the intensity of BBSAN is pro- 

88 portional to the fourth power of /3. The off-design parameter, f3, was defined 

89 by Harper-Bourne and Fisher for convergent nozzles and is extended to the 

90 general case as, 

13 = y'lA/J - Mji ( 1 ) 

91 where is the nozzle design Mach number, which is dependent on the 

92 ratio of the exit area to the throat area, and Mj is the fully expanded Mach 

93 number, which is only dependent on the NPR and the ratio of specific heats 

94 7. It was shown this trend holds over a wide range of fully expanded Mach 

95 numbers for a wide range of convergent and convergent- divergent nozzles. 

96 Note that the 4 th power of (3 may vary slightly due to small effects of heating 

97 and sound emission angle (see Viswanathan et al. [ 16 ] for details). At higher 

98 (3 the relationship is no longer linear and the slope of /3 4 versus NPR drops 

99 off slightly. Equation 1 is relatively independent of TT R. 

100 Harper-Bourne and Fisher [ 10 ] write, “... the intensity of shock noise 

101 is a function only of pressure ratio, and is independent of jet stagnation 

102 temperature and hence jet efflux velocity.” This statement is in the context 

103 of a larger study and conflicts with more recent experimental observations. In 

104 the experiments of Viswanathan [ 16 ] and Kuo et al. [ 7 ] there are noticeable 

105 differences in BBSAN intensity when holding NPR constant and varying 

106 TTR. These differences are often unnoticeable if the jet is non-screeching 

107 (when the jet is heated) compared to a jet that is screeching (when the jet is 

108 cold), which is often the case in small laboratory experiments. This is due to 
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Figure 1: Lossless sound pressure level per unit Strouhal number at R/D = 100 resulting 
from a Md = 1.00, Mj = 1.50, and TTR = 1.00 jet. ^ is the angle from the upstream jet 
axis to the observer centered about the nozzle exit. S 0 is the screech over-pressure. 
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Figure 2: Lossless spectra for a Mj = 1.71 and D = 0.06223 m jet at R/D = 97.86 and 
ip = 90 degrees. The spectra corresponds to a TTR of 1.00, 1.80, 2.20, 2.70, and 3.20. 
This figure is reproduced from Viswanathan et al. [16] with permission. 
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the very large effect that screech tones have on BBSAN (see Andre et al. [ 17 ] 
no for details). It is important to isolate the scaling of BBSAN intensity from 
in the varying NPR and from the effect of various screech intensities, which 
n2 are a function of NPR and TTR. 

n 3 Tam [ 18 ] developed a method for BBSAN prediction and the basic phys- 
ic ical model is described by Tam and Tanna [ 19 ]. Tam argued that the shock 
ns cell structure in the jet could be modeled, following the work of Pack [ 20 ]. 
no The large-scale turbulence in the jet shear layer is modeled as a random su- 
n 7 perposition of instability waves supported by the jet mean flow, as described 
ns by Tam and Chen [ 21 ], Tam [ 22 ] modified the model by Tam [ 18 ] to include 
ns the capability to predict BBSAN from heated jets up to a moderate off-design 

120 parameter, f 3 . A temperature correction factor, T c f, was included to correct 

121 for the over-prediction at all frequencies due to increasing TTR. Tam used 

122 the following empirical correction factor for heated jet predictions, 

Tc/ = fe( 1+1 i^ M 0" 1 (2) 

123 where pj is the fully expanded density, p^ is the ambient static density, and 

124 7 is the ratio of specific heats. Morris and Miller [6] formed an acoustic 

125 analogy for BBSAN and later showed its application to a wide variety of 

126 fully expanded Mach numbers and temperature ratios, for cylindrical, dual- 

127 stream, and rectangular nozzles, with over- and under-expanded jet plumes. 

128 To account for the slight heating effects on the predicted BBSAN relative to 

129 experimental data, Eqn. 2 was used to scale the spectral density. 

130 Recently, Kuo et al. [ 7 ] performed experiments that examined the effects 

131 of heating on BBSAN in the far-held by examining three nozzle geometries. 

132 The first was convergent and the others were convergent- divergent at Md = 

133 1.50 and Md = 1 . 76 . Heating of the jet how was accomplished by simulating 

134 a heated how with a helium-air mixture. Doty and McLaughlin [ 23 ] had 

135 shown that helium-air jets and heated jets have similar physical and acoustic 

136 properties in the far-held. Kuo et al. [ 7 ] examined heating effects for the 

137 Md = 1.50 nozzle at Mj = 1 . 2 , 1 . 4 , 1 . 7 , and 1.9 by varying TTR from 

138 1.0 to 2 . 2 . In the following sections, a scaling relationship is developed for 

139 BBSAN intensity with increasing TTR , and the relationship is compared 
mo with measurement for four cases. The arguments of the scaling relationship 
mi are based upon steady RANS solutions and details of the calculations are 
142 shown. 
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2. Mathematical Analysis 

The Euler equations are rearranged into a linear left hand side operator 
of the LEE and right hand side equivalent sources. The equivalent source 
of the continuity equation is the dilatation and the equivalent source of the 
momentum equation is the unsteady force per unit volume involving veloc- 
ity fluctuations of the shocks and turbulence in the jet plume. The latter 
is of interest for BBSAN prediction. The acoustic pressure is found from 
the convolution integral of the vector Green’s function with the equivalent 
sources. The spectral density is then formulated by the Fourier transform 
of the autocorrelation involving acoustic pressure. The full details of this 
approach are shown in Morris and Miller [ 6 ] and result in, 


S(x,u>) — ploC^ / ... / '%2'%2n™{ x ,y,u)n™( x ,y + rj,u) 

0 ° J —oo n= i m=l W 

xR v nm {y , 77 , r) exp[— ic JT]dTdrjdy 

where S is the spectral density, 7 r” is the n th component of the vector Green’s 
function of the LEE, R^ nn {y ) i 7 , t) is the two-point cross-correlation of the 
equivalent source, x is a vector from the nozzle exit to the observer, and y is 
a vector from the nozzle exit to a source in the jet. rj = r](£, rj, Q is a vector 
between two spatial locations in the jet source region. 

The vector Green’s function of the LEE as shown in Eqn. 3 is defined by 
the solution of, 


and, 


PoKg 9v V gt 
Dt dxi 


S{ x - y)S(t - t)8 0n 


( 4 ) 


D 0 u, 


9* 


Dt 


+ 11,11 a'r, +c2 ^~ = 6<yX ~ y ~~ r)4 


( 5 ) 


where D a is the material derivative about the meanflow and u is the velocity. 
The vector Green’s function is periodic and has the identity it™* ( x , y , lu) = 
Kg{x,y, —u). General analytic solutions of Eqns. 4 and 5 are unknown, 
however, numerical solutions can be found that are related to Lillcy’s [24] 
equation. Strategies to find highly accurate numerical solutions of the vector 
Green’s functions are discussed in Tam and Auriault [25], Raizada [26], and 
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Khavaran et al. [27]. Propagation effects have been examined for BBS AN 
using these techniques by Miller and Morris (2013 IJA) and a ray method 
by Henry et al. [28] . The approach of Miller and Morris is employed here to 
find TTg. R v nm takes the form, 


Km(y ^ T ) = fn(y^)fm(y + 'n^ t + r ) ( 6 ) 

where /? is the equivalent source involving second order fluctuations of the 
momentum term in the governing equations, which is defined as, 


ry d)llfj fh/s?' 

fi= - u ’ 1 dT j - u “s^ 

where u are velocity fluctuations associated with the shocks, s, and the turbu- 
lence, t, and Xi are independent spatial coordinates. In Morris and Miller [6] 
the equivalent source is formulated based on dimensional and physical argu- 
ments involving the speed of sound, c, the integral turbulent length scale in 
the streamwise direction, /, the pressure due to the shock waves, p s , and the 
density, p. In this work we assume that the density and streamwise velocity 
are local instead of ambient values. A model for f%(y,t)fm(y + r/,t + r) is 
formed, 


fn(y^)fm(y + r h t + T ) 


Ps{y)p s (y + y) 

p 2 c 2 l 2 


R{y,y,r) 


( 8 ) 


where R(y, rj , r) is the two-point cross-correlation of the velocity fluctua- 
tions. Assume that the time and spatial terms of R(rj, r) are separable as 
Ribner [29] postulated and model the two point cross-correlation of the fluc- 
tuating turbulent velocity as, 


R(v,V,t) = (l mn K[y) exp [-r 2 /r, 2 ] exp [-(£-« c r) 2 /( 2 ] exp [-(rj 2 + C 2 )/*l] 

(9) 

where a mn are coefficients that can be set for anisotriopic turbulence and K 
is the turbulent kinetic energy. 

Substituting Eqns. 8 and 9 into Eqn. 3 and isolating the integral involving 
r yields, 


exp [— icnr] exp — r 2 /r s 2 exp — (£ — u c t) 2 /1 2 dr. 


( 10 ) 


9 



190 


Integration of expression 10 is performed analytically, 


7 r 


1/2 


exp 


-4g 2 +4iT 2 u e gij--/ 2 T 2 io 2 

4(Z 2 +r|u 2 ) 


^!/r s 2 + M 2//2 


191 Expression 11 is used with Eqn. 3 , 


( 11 ) 


s(x,u) = a 44 , 


/ oo /»oo 3 3 

00 •l'°°n=lm=l 


V 2 , 


X 


a mn /l (y )p s (y )p s (y + r /) 71 " exp 


— 4£ 2 +4ir 2 u c £uj — lg Tg uj 2 
4 Oi+Tfu 2 ) 


( 12 ) 


p 2 c 2 / 2 


^ 1 /r 2 + « 2 // 2 
x exp [-(y 2 + C 2 )//l] dydy. 


192 Over the distance where the spatial correlation is significant we assume 

193 that, 


n™(x, y + rj,u>) — n™(x, y, uj) exp 


. UJ X 

i — i — r ' V 

Cno X 


( 13 ) 


194 as shown by Tam and Auriault [ 25 ]. We now examine the term p s (y + rj) 

195 shown in Eq. 12 . Morris and Miller [6] noted that the variation of the Fourier 
we transform of the shock pressure can be written as, 


/ OO 

Ps{y) exp[i/e,y 1 ]dyi ( 14 ) 

-OO 

197 where k is the spatial wavenumber. It is observed that the variation of 
we p s (ki,r),C) changes little across the jet core and shear layer where the BB- 
199 SAN source is located and is certainly a valid approximation as long as the 
200 variation is small within regions of slowly varying shock pressure. Likewise, 
201 the same argument applies in the untransformed domain in conjunction with 
202 the observation that the spreading rate of the jet is small and that the shock 
203 cell interactions generally occur at the same radius. With these assumptions 
204 it is argued, 


205 


Ps(v + v) -pJv + 0 - 


( 15 ) 


We choose to use Eqn. 15 as it makes the analysis for scaling much simpler. 
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Substituting Eqns. 13 and 15 into Eqn. 12 and isolating the terms of q 
207 and C yields an expression of the integrals involving rj and (, 


f* oo /»oo 


exp 


' — oo J oo 


— 10 JX 2 T] 
Cr^\x\ 


exp 


~ivx 3 ( 


exp 


Hv 2 + C 2 ) 2 


i 2 


( 16 ) 


208 The integrals are evaluated analytically, 


i r/]_ exp 


4 c 2 lad 


( 17 ) 


209 Expression 17 is now used to simplify Eqn. 12 . Let us now restrict our 

210 model to the sideline direction, 6 = it/ 2 . At the sideline direction it 2 is 

211 dominant relative to the other components. 


Se=n/ 2 (u) = 7T 7 Pood 


n 2 *(x,y,u)it 2 (x,y,uj) 


1 — OO J — OO 


exp 


Coo 


P-^y + ^a 22 K( v) — 


-4g L ^4i T;u c &-l-TjJ- 

4(i2+r. 2 «2) 


p 2 c 2 / 2 


Wa 2 + “PI 

~-/ 2 oA~ 


( 18 ) 


x 1 / exp 


4 c 2 

^ L oo J 


dt/dy. 


212 The sources of BBSAN are at relatively discrete locations which is unlike 

213 mixing noise. The integrals of Eqn. 18 are replaced with summations of the 

214 integrand over the source regions, 


A A 


So=*/ 2M = ?r 3/2 /4c4j]^7rf (x,y,uj)n 2 g {x,y, 




x exp 


a = 1 6=1 

i^i^ l Ps(y)p s (y + Q 

Cool X | p 2 C 2 / 2 


a 22 K(y)- 


exp 


-4g 2 -4irfu c gai-ZgT-^ 2 

4(J2+ t 2 u 2) 


^i/r 2 + u 2 // 2 
--l 2 u 2 ' 


x exp 


4 c 2 




( 19 ) 


215 where V) 0 is the local source volume and distance £ around each shock wave 

216 shear layer interaction, a and b. The total number of shock wave shear layer 
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interactions is A. a 22 is an element of the a mn tensor. If we restrict onr 
analysis to the contribution from a single shock wave shear layer interaction 
then £ = 0 and Eqn. 19 becomes, 


Se=ir/2(u) — 


r 3/2 2 

r OO 

prefactor propagation 


* 7T“-7 Tg exp 





1 

— to 

£ 

to 

exp 

4c 2 

L ^ 00 J 


Ws + T 2 U 2 ) 

a 22 p 2 s l\K 


x 


P 2 C 2 V 2 ^1/t 2 + 


V 


( 20 ) 


where V is the local source volume. The hrst term, ‘prefactor, ’ has no effect 
on the scaling of BBSAN while varying TTR. The second term, ‘propaga- 
tion,’ is the vector Green’s function components. It quantifies the effect of 
the sound propagation of BBSAN and is important for capturing tempera- 
ture effects. The third labeled term, ‘source,’ results from the choice of the 
equivalent source. The turbulent kinetic energy, K, and the local properties 
of p, c, and the integral scales of turbulence at the shock wave shear layer in- 
teraction control the scaling in the ‘source’ term. The two exponential terms 
of Eqn. 20, based on numerical variation relative to the ‘source’ term, have 
little effect on the variation of the spectral density with increasing tempera- 
ture. However, the exponential terms are included in the predictions in the 
following section for accuracy and completeness. 

Morris and Boluriaan [30] have shown that \7r 2 (x, y\ uj) \ 1 = uj 2 / (lO^c^o; 2 ) 
for the far-held at the sideline angle for axisymmetric jets. Here it results in, 


Sn= 


0=7r/2 


,CJ = 


p 2 00 a 2 2Kp 2 uj 2 


x exp 


I 67 T 1 / 2 cl c x 2 p 2 c 2 ^/1 /t 2 + u 2 /l 2 s 
~ —1 2 t 2 u 2 


Wl + 7>c) 



1 

— to 

£ 

to 

exp 

4c 2 

L ^00 J 


V 


( 21 ) 


Equations 20 and 21 yield the BBSAN intensity from a single shock wave 
shear layer interaction in the sideline direction. In the following sections we 
will evaluate Equation 19 for the acoustic intensity from BBSAN at the peak 
frequencies in the sideline direction for multiple jet conditions for all source 
interactions. 
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M d 

D (m) 

Mj 

NPR 

Origin 

Type 

1.00 

0.0508 

1.50 

2.42 

Bridges and Brown 

Convergent 

1.50 

0.0127 

1.20 

3.67 

Kuo et al. 

Convergent- D ivergent 

1.50 

0.0127 

1.70 

4.94 

Kuo et al. 

Convergent- D ivergent 

1.50 

0.0127 

1.90 

6.70 

Kuo et al. 

Convergent- D ivergent 


Table 1: Properties of the jet flows. For each row a simulation is performed at TTR 1.00 
to 2.50 at increments of 0.10 and from TTR 2.50 to 4.00 in increments of 0.25. In total 
80 steady RANS solutions are performed. 


239 3 . Results 

240 Four cases are selected to exercise Eqn. 19. The cases represent over- 

241 expanded and under-expanded conditions over a range of Mach numbers for 

242 a convergent and convergent-divergent nozzle. The nozzles and operating 

243 conditions are shown in table 1 and the TTR varies from 1.00 to 4.00 for 

244 each case. The first row of the table shows the conditions of the convergent 

245 nozzle and has corresponding data collected from the SHJAR experiment of 

246 Bridges and Brown [ 8 ]. The remaining three rows of the table correspond to 

247 three of the four conditions performed in the experiment of Kuo et al. [7]. 

248 A CFD calculation is performed for each experimental condition summa- 

249 rized in table 1. The arguments of the acoustic analogy are related to the 

250 steady RANS solution. The equivalent sources could easily be informed by 

251 a more advanced simulation that uses LES or simpler empirical models. 

252 3.1. Steady Reynolds- Averaged Navier- Stokes Solutions 

253 The Wind-US CFD (see Nelson [31] for details) solver is used to calculate 

254 the steady RANS solutions. Calculations are performed from TTR = 1.00 

255 to 4.00 in increments of 0.1 for TTR = 1.00 to 2.50 and increments of 0.25 

256 from TTR = 2.50 to 4.00. All simulations are axisymmetric and are closed 

257 by the Menter [32] Shear Stress Transport (SST) turbulence model. Details 

258 for these types of simulations and experimental validation of the flow-fields 

259 have been discussed by Miller and Veltin [33]. 

260 By examining table 1, it can be shown that only four NPR are required 

261 and only four unique values of (3 result. Changing TTR while holding the 

262 NPR constant results in nearly identical shock-cell structures. Contours of 

263 p/Poo are shown in Fig. 3 where the axes are normalized by the jet diameter. 
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264 Steady RANS solutions are mirrored about the x-axis for illustration pur- 

265 poses. The jet conditions in Fig. 3 are a) M d = 1.50, Mj = 1.20, TTR = 1.00, 

266 D = 0.0127 m, b) M d = 1.00, Mj = 1.50, TTR = 1.00, D = 0.0508 m, c) 

267 M d = 1.50, Mj = 1.70, TTR = 1.00, D = 0.0127 m, and d) M d = 1.50, 

268 Mj = 1.90, D = 0.0127 m, TTR = 1.00. The circles in parts a) through 

269 d) represent the time averaged locations where conical oblique shock waves 

270 interact with the jet shear layer. The shock cell shear layer interactions rep- 

271 resent the positions where BBSAN sources are located. At each shock wave 

272 shear layer interaction, the held variables are extracted as a function of TT R 

273 from the steady RANS solutions. 

274 To illustrate the relative source strength and location of the BBSAN 

275 sources, a numerical investigation is performed with a M d = 1.00, Mj = 1.50, 

276 TTR = 1.00, and D = 0.0254 m jet. A steady RANS solution of this how- 

277 held is shown in Fig. 4. Part a) shows contours of shock pressure and part 

278 b) shows contours of turbulent kinetic energy. Part c) shows contours of the 

279 integrand of the model of Morris and Miller [6] at ^ = 90 degrees and R/D = 

280 100. The contours of part c) represent the relative strength of BBSAN at the 

281 peak BBSAN frequency. At this source location the high values of p s and K 

282 can be observed in part a) and b) respectively and indeed, as theory suggests, 

283 correspond to shock wave shear layer interactions. Measurements, such as 

284 those of Norum and Seiner [34], show that the BBSAN source is further 

285 downstream than prediction, however, in their study there is no account for 

286 the refraction effects of the jet shear layer. In supersonic jets, refraction 

287 effects can make sources appear multiple diameters downstream from their 

288 actual location. 

289 3.2. Scaling of Broadband Shock-Associated Noise with Temperature 

290 A comparison of the predicted source scaling of the BBSAN with mea- 

291 surement is shown in Fig. 5 for the convergent nozzle, Mj = 1.50, and varying 

292 TTR. The contribution predicted by Eqn. 19 is shown as a black line with 

293 round circles. The dashed line with triangles is the prediction of Morris and 

294 Miller [6] without T c f. The data from the Small Hot Jet Acoustic Rig (SH- 

295 JAR) experiment is shown as red squares. Experimental values represent the 

296 maximum BBSAN at the sideline location of the jet. The evaluation of the 

297 intensity has been performed based on the locations shown in Fig. 3 part 

298 b). The factor T c f is not included in the predictions using Eqn. 19 or those 

299 derived from Morris and Miller [6]. It is clear that with increasing TTR the 
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Figure 3: Contours of pressure p/poo of the four jet families studied. The circles represent 
locations of the oblique shock wave shear layer interactions. Flow-field data is extracted 
as a function of TTR at these locations. The jet conditions shown are, a) M d = 1.50, 
Mj = 1.20, TTR = 1.00, b) M d = 1.00, M 3 = 1.50, TTR = 1.00, c) M d = 1.50, M 3 = 1.70, 
TTR = 1.00, and cl) M d = 1.50, M 3 = 1.90, TTR = 1.00. 
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Figure 4: A Md = 1.00, M. ; = 1.50, TTR = 1.00, D = 0.0254 m jet produces contours of, 
a) shock pressure, b) turbulent kinetic energy, K, c) integrand of the model of Morris and 
Miller [6] at the sideline location at R/D = 100 and peak BBSAN frequency (contours of 
BBS AN source strength). Note in c), the maximum BBSAN source occurs at the location 
where the shock wave shear layer interactions occur. 
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300 prediction of Eqn. 19 initially increases linearly with TTR and eventually 

301 saturates. 

302 The following comparisons are performed with the experiment of Kuo et 

303 al. [7] who used a M d = 1.50, D = 0.0127 m nozzle, and varied the TTR. 

304 Figure 6 shows comparisons at R/D = 100 and at the sideline location 

305 with Mj = 1.20. These particular jet conditions produced no screech tones 

306 through their entire temperature range. The scaling of the source shows a 

307 clear increase from TTR = 1.00 and saturates at relatively the same rate as 

308 the experiment. 

309 For the next comparison, the same nozzle and observer location is retained 

310 but the jet operates at Mj = 1.70. Comparisons between the predicted peak 
3n BBSAN and Kuo et al. [7] are shown in Fig. 7. Unlike the previous case, these 

312 jet conditions produced very strong screech tones. The over-pressure of the 

313 screech, S 0 . is marked at each data point in the figure. S Q is a measure of the 

314 maximum screech amplitude minus the broadband level at the fundamental 

315 screech frequency. The screech frequency is often lower than the peak BBSAN 

316 frequency. To illustrate this point, reexamine Fig. 1 at the sideline location, 

317 where the fundamental screech tone frequency is lower than the peak BBSAN 

318 frequencies. The tone labeled screech is the fundamental screech tone and 

319 its overpressure is approximately 12 dB. In Fig. 7 the maximum BBSAN 

320 has corresponding values of S a that change with TTR. At low temperatures 

321 the screech over-pressure is large and as the TT R increases the screech over- 

322 pressure approaches zero. In Fig. 5 where S 0 is relatively constant and non- 

323 zero and Fig. 6 where there is no screech, the agreement between prediction 

324 and experiment is arguably better. If the screech over-pressure is constant 

325 through the range of TTR or preferably, zero, as it is typically in full-scale 

326 engines unlike small nozzles, then the effects of screech on the BBSAN are 

327 relatively the same and the scaling of BBSAN is accurately captured. In 

328 Fig. 7, one may observe the correct trend of BBSAN saturation starting at 

329 TTR = 1.4. The BBSAN is amplified a large amount due to the large screech 

330 amplitudes present at low temperatures. 

331 A final comparison is made between the prediction of the scaling using 

332 Eqn. 19 and the experiment of Kuo et al. [7] in Fig. 8. The fully expanded 

333 Mach number is increased to 1.90 and all other conditions are retained. The 

334 similarity between the prediction of Eqn. 19 and measurement of Kuo et 

335 al. [7] at moderate to higher TTR is similar to that of Fig. 7. At low temper- 

336 atures the screech tones have disrupted the trend due to reasons previously 

337 described, and the correct scaling of BBSAN is not captured. If screech were 
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Figure 5: Variation of maximum BBSAN intensity from a M ( i = 1, Mj = 1.5 (experiment 
performed at Mj = 1.469) jet relative to increasing TTR. The observer is located at 
R/D = 100 and ip = 90 degrees. The prediction of Eq. 19 is compared with measurement 
of the maximum BBSAN. 
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Figure 6: Variation of maximum BBSAN intensity from a Aid = 1.5, Alj = 1.2 jet relative 
to increasing TTR. The observer is located at R/D = 100 and if> = 90 degrees. The 
prediction of Eq. 19 is compared with measurement of the maximum BBSAN. 
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Figure 7: Variation of maximum BBSAN intensity from a Aid = 1.5, Alj = 1.7 jet relative 
to increasing TTR. The observer is located at R/D = 100 and if> = 90 degrees. The 
prediction of Eq. 19 is compared with measurement of the maximum BBSAN. 
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Figure 8: Variation of maximum BBSAN intensity from a Md = 1.5, Mj = 1.9 jet relative 
to increasing TTR. The observer is located at R/D = 100 and i/j = 90 degrees. The 
prediction of Eq. 19 is compared with measurement of the maximum BBSAN. 


not present within the experiments then the trends at low TTR for Figs. 7 
and 8 will yield a lowered BBSAN amplitude, and eventual saturation as 
seen in Figs. 5 and 6 will occur. 

It is evident that the inclusion of using local properties for the streamwise 
velocity component and density (in the denominator of the source term), 
instead of ambient quantities, and the combination of the vector Green’s 
function that is amplified by the shear layer, yields a model that is more 
consistent with experiment. 

Unheated and slightly heated jets are difficult to predict due to the rapid 
variation of BBSAN intensity. Not only are the predictions very difficult 
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to conduct both mathematically and in terms of implementation but the 

349 experiments are very difficult to perform, especially with nozzles on the order 

350 of 10 -2 meters. The difficulty of acquiring excellent experimental data cannot 

351 be overstated. 

352 Screech tones are extremely sensitive to laboratory conditions and are 

353 highly nonlinear, however, all fluid dynamic phenomena are deterministic. 

354 When the screech over-pressure is very high at low TTR , the BBSAN is 

355 lowered in frequency and raised in amplitude. It could be possible for an 

356 experiment to be performed for the same jet conditions as shown in this 

357 paper, but without screech and without screech’s effect on the mixing noise 

358 or BBSAN. It is expected that the peak BBSAN intensity levels, without the 

359 influence of the discrete tone, will compare with the developed theory. 

360 4. Conclusion 

361 BBSAN intensity saturates with increasing jet stagnation temperature. 

362 This saturation occurs due to the balance between the source term and the 

363 propagation effects. An equivalent source for BBSAN is proposed that takes 

364 into account the scaling of both NPR and TTR. The scaling term is con- 

365 tained within the developed acoustic analogy and contains the effects of the 

366 equivalent source and propagation separately. This acoustic analogy is eval- 

367 uated with arguments corresponding to four families of disintegrated jets. 

368 Evaluation involves extracting the local properties at the shock wave shear 

369 layer interactions from steady RANS solutions. These local field variables are 

370 arguments of the source term. Comparisons of the predicted peak BBSAN 

371 intensity from shock wave shear layer interactions show the same trend as 

372 measurement. The predictions, like the experiments, show eventual satura- 

373 tion with increasing jet stagnation temperature. At very high temperature 

374 ratios predictions show that saturation ceases and the BBSAN intensity will 

375 again rise. Higher fidelity measurements, without screech, are required to 

376 further validate this theory. 
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